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Quantum chromodynamics is notoriously difficult to solve at nonzero baryon density, and most 
models or effective theories of dense quark or nuclear matter are restricted to a particular density 
regime and/or a particular form of matter. Here we study dense (and mostly cold) matter within 
the holographic Sakai-Sugimoto model, aiming at a strong-coupling framework in the wide density 
range between nuclear saturation density and ultra-high quark matter densities. The model contains 
only three parameters, and we ask whether it fulfills two basic requirements of real-world cold and 
dense matter, a first-order onset of nuclear matter and a chiral phase transition at high density to 
quark matter. Such a model would be extremely useful for astrophysical applications because it 
would provide a single equation of state for all densities relevant in a compact star. Our calcula¬ 
tions are based on two approximations for baryonic matter, firstly an instanton gas and secondly 
a homogeneous ansatz for the non-abelian gauge fields on the flavor branes of the model. While 
the instanton gas shows chiral restoration at high densities but an unrealistic second-order baryon 
onset, the homogeneous ansatz behaves exactly the other way around. Our study thus provides all 
ingredients that are necessary for a more realistic model and allows for systematic improvements of 
the applied approximations. 


I. INTRODUCTION 
A. Motivation 

Quantum chromodynamics (QCD) at nonzero baryon densities and zero (or very small) temperatures is very difficult, 
and progress has mostly been made with the help of various effective theories and phenomenological models. What we 
know for certain is that there is a first-order transition from the vacuum to nuclear matter at nuclear saturation density 
n o ~ 0.16 fm -3 and that QCD becomes weakly interacting at ultra-high densities due to asymptotic freedom HH2], 
wherefore sufficiently dense matter is a weakly coupled gas of quarks. Hence there must be chiral and deconfinement 
phase transitions presumably at moderate densities well in the strong-coupling regime - from nuclear to quark 
matter 1 . All current rigorous studies are restricted to a particular density regime, for example, “ordinary” nuclear 
physics at and slightly above saturation density, lattice studies for baryon chemical potentials smaller than or at 
most of the order of the temperature, and perturbative studies at extremely large densities. But even giving up on 
some rigor, it is very difficult to come up with a reasonable model that describes QCD matter over a wide density 
regime. For example, the quark-meson model, see for instance Refs. [3 Ho], and the Nambu-Jona Lasinio (NJL) model 
HUE], see for instance Refs. |l3til5j . and variants thereof can be very useful to get some insight into quark matter 
phases and possibly the chiral and deconfinement phase transitions, but they usually do not include nuclear matter, 
although it is possible HU. On the other hand, nucleon-meson models, such as in Refs. HM5, are based on the 
properties of nuclear matter at saturation and may describe moderately dense nuclear matter realistically, but give 
a poor description, if at all, of chirally restored quark matter. Therefore, attempts to study strongly coupled dense 
matter from the onset of baryons all the way to quark matter at high densities are mostly based on patching together 
various models and, not surprisingly, depend on many parameters [Hi [35]. Having a single model at hand - and be it 
only a very distant relative of the fundamental theory - would be very interesting not only from a theoretical point 
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1 Although most models predict a first-order phase transition for small temperatures, it is conceivable that there is a smooth crossover 
[3H5 ] from nuclear matter to color-flavor locked (CFL) quark matter [21E]- In this paper, we shall ignore color-superconducting phases 
such as the CFL phase, although it is an interesting question whether it can be implemented in the Sakai-Sugimoto model, possibly 
following existing attempts within holography [§]. Also, we shall work in the chiral limit, such that the transition between chirally 
broken baryonic matter and chirally restored quark matter must be a phase transition in the strict sense. 
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of view, but also for the applications in the context of compact stars. Compact stars have a density profile, and thus 
knowledge about the strongly coupled physics from nuclear saturation density possibly up to quark matter densities 
in the center of the star is required [ 21] . 


B. Model 

The gauge-gravity correspondence [25U27] provides a powerful tool to study strongly coupled physics. Although 
the existence of a gravity dual of QCD is suggested by general principles, it is currently unknown and probably out of 
reach in the near future, and most studies have focused on theories that are different, sometimes very different, from 
QCD. The most prominent example is Af = 4 super Yang-Mills theory where Nf quark flavors are introduced in the 
gravity dual by Nf D7-branes in the background given by N c D3-branes; for studies in this setup aiming at the QCD 
phase diagram see for instance Ref. [28] and, including baryonic matter, Refs. [251120] ■ The Sakai-Sugimoto model 
[3RMI is the holographic model that currently comes closest to QCD. In this paper we ask the question whether 
it can be employed in the dense matter context we have just outlined. Like the D3/D7 system, the Sakai-Sugimoto 
model is a “top-down” approach, i.e., it is rigorously based on an underlying type-IIA string theory, and the dual 
field theory is known. Namely, in a certain limit, it is dual to large- N c QCD. Here, we apply several approximations 
and extrapolations, giving up some of this rigor in favor of feasibility and of potentially coming closer to real-world 
N c = 3 QCD at large baryon chemical potentials: firstly, as in most previous applications of the model, we work in 
the classical gravity approximation, which is a good approximation for very large values of the ’t Hooft coupling A 
(although we shall extrapolate some of our results to small A), while large- N c QCD is the dual field theory only in the 
opposite, inaccessible, limit of small A. Secondly, our main results are obtained in the “deconfined geometry”, which 
has a richer phase structure concerning chiral symmetry. While the “confined geometry” can be connected to the 
confined phase of the dual field theory, this is less clear for the deconfined geometry . We thus do not a priori 

know the (phases of) the field theory we are actually working in, especially regarding confinement/deconfinement 
(chiral symmetry, on the other hand, and its spontaneous breaking are implemented unambiguously). Thirdly, as 
in most related studies, we shall work in the probe brane limit, i.e., we shall not include the backreaction of the 
flavor branes on the background geometry (for studies going beyond this quenched approximation, see Refs. [57] 1581). 
And, we shall make use of a certain, not uniquely defined, prescription for the non-abelian Dirac-Born-Infeld (DBI) 
action to all orders in the string tension. Therefore, at very large densities, where effects of backreaction as well as 
higher-order terms in the string tension may become important, our results must be considered as an extrapolation 
from the low-density regime, where the approach is rigorous. 


C. Goal and relation to previous works 

One important benefit of the given holographic approach is that there is a well defined way to implement baryons. 
Following the general idea of a “baryon vertex” in AdS/CFT [39] 00], baryons in the Sakai-Sugimoto model are 
implemented as D4-branes, which have N c string endpoints, wrapped on the 4-sphere of the background geometry. 
Equivalently, they can be considered as instantons of the gauge field theory on the connected D8- and D8-branes 
[5710T1 - 05| . Baryon properties have been computed with the help of the Belavin-Polyakov-Schwarz-Tyupkin (BPST) 
instanton solution [441 . using instantons in flat space as a trial function for the equations of motion of the curved 
geometry [42 ], although improvements to this approach might be necessary [45H45] , Here we are not interested in 
single baryons in vacuum, but in baryonic matter. We shall restrict ourselves to homogeneous matter; for studies of 
inhomogeneities and instanton crystals in the confined geometry of the Sakai-Sugimoto model, see Refs. [501152] . A 
simple and very useful approximation for homogeneous baryonic matter makes use of pointlike instantons [53s and has 
been employed in both confined and deconfined geometries. This approximation shows an unrealistic second-order 
phase transition from vacuum to baryonic matter (in both geometries), and, also unrealistically, low-temperature 
baryonic matter remains the preferred phase for arbitrarily large chemical potentials. In the confined geometry, this 
is obvious from the geometry of the model in the probe brane approximation and a consequence of the \&rge-N c 
limit. (Although it might be possible to see signatures of chiral restoration also in the confined geometry [53], in 
particular by going beyond the probe brane approximation [55] .) In the deconfined geometry, where chiral restoration 
at large chemical potentials is a priori possible, it turns out that the free energy of the baryonic phase asymptotically 
approaches that of the chirally restored phase, but always remains lower [55j . The main question we address in this 
paper is whether going beyond the pointlike approximation of baryons in the deconfined geometry can remedy one or 
both of these unphysical properties. It has been argued that, at least in the confined geometry, finite-size instantons 
indeed give rise to a first-order baryon onset [S2 EZj- The first part of our calculations is based on this finite-size 
instanton approach, but we shall argue that a fully dynamical calculation does not show a first-order onset, see Sec. 
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III D and in particular footnote [3J The second part of our calculations makes use of a simple homogeneous ansatz 
instead of the instanton solution, which, for the confined geometry, has been introduced in Ref. |58) . 


D. Outline of the paper 


We start with presenting the general setup in Sec. [TTJ where we introduce the gauge field action on the flavor 
D8-branes. The next two sections deal with our two different approaches: in Sec. m we discuss the instanton gas 
and Sec. |IV| is devoted to the homogeneous ansatz. In the case of the instanton gas, we discuss the expansion of the 
Dirac-Born-Infeld (DBI) action for small non-abelian field strengths and the relation to the pointlike approximation in 
Sec. |III A| Our actual calculation is then performed using the full DBI action, and we derive the necessary equations 
and explain the numerical evaluation in Secs. |III B| and |III C| These sections are rather technical and readers only 
interested in the result may jump to Sec. HID where we present our results (including the results for the confined 
geometry). Sec. IV is structured similarly: we first explain our approach and the calculation in Sec. IV A[ before we 
present the numerical results in Sec. |IV B] The main results are Figs.[3 and [4] for the instanton gas and Figs. [6] and [7] 
for the homogeneous ansatz. Finally, we give our conclusions in Sec. [V 


II. SETUP 

We start by setting the stage without going into details of or reviewing the Sakai-Sugimoto model since there are 
numerous works in the literature where this has been done, see for instance Refs. (55. 5Bi, 5?11551 - 1521 . In the following 
we focus on the features of the model that are relevant for us and write down the action we are using. 

The Sakai-Sugimoto model incorporates a deconfinement phase transition by allowing for two different background 
geometries, and realizes a chiral phase transition by two different embeddings of the flavor branes in these background 
geometries. In the original version of the model, the flavor D8- and D8-branes, corresponding to left- and right- 
handed fermions, are maximally separated at the boundary U = oo, where U is the holographic coordinate. Denoting 
this asymptotic separation by L , maximal separation means L = tt/Mkk, where the Kaluza-Klein mass Mkk is 
the inverse radius of the compactified extra dimension X4 of the model, X4 = X4 + 2 tt/Mkk- In this version, 
the chiral and deconfinement phase transitions coincide, i.e., in the confined geometry the D8- and D8-branes are 
always connected, which corresponds to the chirally broken phase, and in the deconfined geometry they are always 
disconnected, indicating an unbroken chiral symmetry. Both phase transitions occur at a critical temperature T c = 
Mkk/(27t) and do not depend on the baryon chemical potential [i if the backreaction of the flavor branes on the 
background geometry is neglected. The reason is that the chemical potential is introduced via the gauge fields on 
the flavor branes and thus - within this probe brane approximation - cannot affect the background geometry and 
therefore has no effect on either deconfinement or chiral phase transitions. 

In this paper, we are mostly interested in non-antipodal asymptotic separations where L <C tt/Mkk- This version 
of the model can be considered as the “decompactified” limit, because it is achieved by either choosing L very small 
compared to a fixed radius MjT^, or by keeping L fixed and choosing a very large radius, i.e., by decompactifying the 
extra dimension. Of course, the dual field theory is only truly four-dimensional in the opposite limit, i.e., for a very 
small radius of the extra dimension, and effects of the fifth dimension can be expected to become important in the 
decompactified limit. Since a large radius corresponds to a very small Kaluza-Klein mass, the critical temperature 
for deconfinement goes to zero in the decompactified limit. Thus, when we work in the deconfined geometry, we shall 
ignore that below T c the confined geometry is preferred. Besides the Kaluza-Klein mass Mkk and the asymptotic 
separation of the flavor branes L , the only other free parameter of the model is the’t Hooft coupling A. 

For our purpose, the decompactified limit is most interesting because here the chiral phase transition is no longer 
locked to the deconfinement phase transition and depends on T and /i. The resulting phase diagram in the T-ji plane 
in the absence of baryons [63j (and including a background magnetic field [641 5511 shows striking similarities to that 
of a Nambu-Jona-Lasinio (NJL) model, in accordance with Ref. ;66]t57], where the connection to the NJL model was 
first pointed out. Two differences to the NJL model are the current quark masses, which are very naturally included 
in NJL but very difficult to implement in Sakai-Sugimoto (although some attempts have been made [581171] ), and the 
order of the chiral phase transition which, in the chiral limit, can be second or first order in NJL but is always first 
order in Sakai-Sugimoto. 

The asymptotic separation of the flavor branes L therefore extrapolates between the original version of the Sakai- 
Sugimoto model, where the dual theory is large- N c QCD (at least in the inaccessible limit of small A), and the 
decompactified limit, which has an NJL-like dual. Since dense matter at large N c is very different from that at N c = 3 
mm, we are less interested in the more rigorous large-1V C limit and focus our attention on the NJL-like limit, even 
though some important physics related to the gluon dynamics are thrown out. One may for instance ask what it really 
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means to include baryons in a limit of the model where there is no confinement. We shall not discuss this question 
in detail and simply observe that instantons and thus baryon number can be introduced in the deconfined geometry 
completely analogously to the confined geometry. 

The two scenarios in which we study baryonic matter are 

• Confined geometry with antipodal separation, L = (free parameters: A, MrkJ- This scenario is used as 

a reference and preparation for the more complicated calculation in the deconfined geometry. The calculation is 
simpler than that of the deconfined geometry because the embedding of the connected branes does not have to 
be determined dynamically. All equations and derivations are deferred to the appendix; in the main part we only 
present the results, which are also useful as a comparison to the existing literature. The results are independent 
of temperature, and thus the only relevant thermodynamic variable is the baryon chemical potential. 

• Deconfined geometry (free parameters: L, A, Mkk)- Our main results are obtained in this geometry and all 
derivations in the main text refer to this case. The results depend on temperature and baryon chemical potential. 
Although we have in mind the decompactified limit when we discuss this geometry for all temperatures, we do 
not make use of L <C 7t/Mkk in the actual calculation and simply treat L as a free parameter. 

Our calculation starts from the gauge field action on the D8- and D8-branes, which consists of a Dirac-Born-Infeld 
(DBI) and a Chern-Simons (CS) contribution, 


S = S'dbi + Scs ■ 


(1) 


We first discuss the DBI term, which has the form 


r 1 / T r r°° _ 

Sdbi = 2T 8 V 4 / dr d 3 X / dU e~* ^detfo + 2 -Kot’T ), 

JO j JlJr 


( 2 ) 


where a! = £ 2 is the string tension with the string length £ s , T§ = l/[(27r) 8 t'®] is the D8-brane tension, and e 4, = 
g s (U/R ) 3 / 4 is the dilaton with the string coupling g s and the curvature radius of the background geometry R. We 
have already performed the trivial integration over the 4-sphere, resulting in the volume factor V 4 = 87r 2 /3. The 
remaining integral is taken over Euclidean time r, 3-dimensional position space, and the holographic coordinate U, 
from the tip of the connected D8- and D8-branes U = U c to the holographic boundary U = 00 . This integral is taken 
over one half of the connected branes, hence the prefactor 2 to account for both halves. The induced metric on the 
flavor branes g in the deconfined geometry is given by 


ds 2 


D8 


= ( D 


3/2 


(f T dr 2 + Ap/.Vb/.V'j + - 


3/2 


{duXi)2 


dU 2 + U 2 dn'l 


(3a) 


where i = 1,2,3, and dfl 2 is the metric of the 4-sphere. The function X^U) describes the embedding of the flavor 
branes in the background geometry, and 


U 3 

f T = l— U t 
lT ~ U 3 ’ 


(4) 


where Ut is the location of the tip of the cigar-shaped t-U subspace, related to temperature and curvature radius via 

(5) 


3 U 1/2 
T = 6 u t 


4t r R 3 / 2 ' 

The following relations between the parameters of the model will be useful, 

2 R 3 Mkk 


R = ng s N c t 


3 

S ’ 


A = 


P. 


( 6 ) 


where N c is the number of colors, i.e., the gauge group of the dual field theory is SU(N C ). 

We are considering two flavors, i.e., a U( 2) gauge theory in the bulk, whose local symmetry group corresponds to 
a global symmetry in the dual field theory. The field strength tensor is decomposed into a U{ 1) and an SU{ 2) part, 


Ffiv — Ffiv + 1 Fpv — F(f y a a , (7) 

where p, v = 0,1, 2,3, U, and the Pauli matrices tr a (a = 1, 2, 3). We shall use the same notation for the gauge fields, 
i.e., and A M for abelian and non-abelian components, respectively. The physics we are interested in requires us 
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to introduce a baryon chemical potential and baryons. The quark chemical potential (related to the baryon chemical 
potential by a factor N c ) is introduced as the boundary value of Ao, while the baryons are introduced in the SU(2) 
part via the field strengths Fjj, Fuj. Thus our only non-vanishing field strengths will be F 0 u = —duA 0 , Fij, Fm- We 
shall set i 7 )o = diA$ = 0 from the beginning because we are only interested in homogeneous matter, i.e., A$ will only 
depend on the holographic coordinate U , not on X. 

Within this ansatz, the CS contribution is 

n r r°° 

Scs = dr J d 3 X J dU A { ;\\[F %jl F ku }e llk . (8) 


III. BARYONIC MATTER FROM AN INSTANTON GAS 


A. Expansion in non-abelian field strengths 


The DBI action as written in Ecp ([2]) does not define how to treat the case of non-abelian gauge fields. A general 
form, keeping all orders of the non-abelian gauge fields, is not known. In expansions for small F (more precisely, 
for small string tension a'), usually the symmetrized trace is used for all terms of 0(E 4 ) and higher, which however 
is known to be incomplete starting from 0(F e ) [741 . For our instanton gas approximation, we discuss two different 
approaches: in the current section, we expand the DBI action for small non-abelian field strengths up to quadratic 
order, while keeping the abelian field strength to all orders, following Ref. [56]. Later, in Sec. IIIB we shall employ 
an “abelianized” prescription that keeps abelian and non-abelian field strengths to all orders and that we use for our 
numerical calculation (for our purposes, it is actually simpler not to expand the square root). 

For the expansion to second order in the non-abelian field strengths we need to replace the square root of the DBI 
action ([2]) as follows, 


a/ det(g + 27 tol'F) —¥ \J det(ff + 2ira'F) |l — ^ (g + 27ra'F) 1 F | + (D(F 4 ) 


(9) 


with 


det (<7 + 2na'F) = U 4 I — 


i ?\ 3/4 


U 


1 + 11 3/tI? - «f 2 


0 ) 


(10a) 


Tr 


(g + 2ira'F)~ 1 F 


-[§) nm- 


2/tTi-[F. 


rtj\ 


1 + u*f T x'Z - Qq 


(10b) 


where the trace on the left-hand side of Eq. (10b I is taken over 8 + 1-dimensional metric space and 2-dimensional 
flavor space, while the trace on the right-hand side is only taken over flavor space. We have replaced Ao —» iAo, which 
is necessary since we work in Euclidean space-time, and the prime denotes derivative with respect to u. Moreover, 
we have introduced the dimensionless quantities it, X 4 , ao, defined in Table [I] together with all other dimensionless 
quantities used in this paper 2 . 

Inserting these results into the DBI action Q, we obtain 


Sdbi = Af J dr J d 3 X J du n 5,/2 \Jl+ vA/tx'I — dg 


1 + Tr[F: 


4 U 3 


h (2W) 2 Tr[F 2 / ] 


2 1 + u 3 /tx 2 — a'o \ 


where we have abbreviated 


A7 = 


2 T 8 V 4 


R 5 (M kk R ) 7 = 


N c R 2 (M kk R ) 7 
67 r 2 (27ra') 3 


( 11 ) 


( 12 ) 


2 There is a slight difference in convention compared to Refs. |53l I55|, which otherwise use the same notation: in that references, for 
example do = 27ra / / RAq. We have included additional powers of the dimensionless factor Mkk^ because then all observable quantities 
(chemical potential, baryon density, temperature) are directly obtained from their dimensionless counterparts by choosing values solely 
of A and Mkk (and not also values for R and a !, which have no direct interpretation in the dual field theory). 
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X 4 ,£ 

U, Uc,UT,UKK,Z,p 

t 

9(2) 

00, p, h 

ni 

Mkk 

1 

1 

R{M kk R) 2 

‘I'KOL 4 - 7 T 

N c R(MkkR) 2 96 tt 4 

R(MkkR) 2 

Mkk 

R(MkkR) 2 AMkk 

X 2na' A 2 M£ k 


TABLE I: Relation between the dimensionless quantities (first row) and their dimensionful counterparts, for example aq = 
MkkA^jI = T/Mkk, etc. Generally, we use capital letters for dimensionful and small letters for dimensionless quantities. 
Exceptions are p (here we use only one symbol for both quantities, which does not cause confusion since both quantities 
appear always well separated from each other), q(z) [where the dimensionful version is denoted by q(Z)], p (for which we never 
introduce a dimensionful version), and n/ (where the dimensionful version is Ni/V, Ni being the instanton number, V the 
3-volume, and ni the dimensionless instanton number density). 


Instead of solving the full equations of motion for abelian and non-abelian fields, we shall for simplicity employ 
an ansatz for the non-abelian field strengths that is based on the BPST instanton solution. We then will have to 
minimize the free energy with respect to the parameters that are introduced by this ansatz (the number density of 
instantons and their width), together with solving the equations of motion for the abelian held fio and the embedding 
function X4. The instanton solution is best introduced in a coordinate that extends continuously over both halves of 
the connected flavor branes. This new coordinate Z is defined via 


u = <u’ + uW 

( c c ’ ’ dZ 3U 1 / 2 


where 


fc = 1- 


ul 

u 3 


(13) 


(14) 


While U £ [U c , 00 ], we have Z £ [— 00 , 00 ], where Z = 0 corresponds to the tip of the connected branes. The in stanton 
solution is well known for the hat-space Yang-Mills (YM) action. As explained in appendix [Aj see Eq. (All I, within 
this solution the traces of the held strengths become 


Tr [F? z ] = 


6 4 (p/7) 4 

7 2 K 2 + 0/7) 


214 > 


Tr[^.] = 12 


4 (p/7) 4 


\e + {ph) 


214 ’ 


Tr [FijF k z\eijk = - 


12 4 (p/7) 4 

7 [£ 2 + (P/7) 2 ] 4 


where we have used Tr[cr a cr a ] = 6 , abbreviated £ 2 = (X — X 0 ) 2 + \{Z — Z 0 )/ 7 ] 2 , and 


3 U, 


3/2 


3 u. 


3/2 


7 = 


2 R 3 / 2 


-(Mkk-R ) 3 


(15) 


(16) 


The instanton is located at the point A ' 0 in position space and the point Z 0 in the bulk and its width is characterized 
by p. The factor 7 was introduced in the BPST solution (A 8 ): while the instanton width in the holographic coordinate 


is p, it is p /7 in position space. In the confined geometry with maximally separated flavor branes, 7 is just a constant, 
see Eq. (AlOl, and one can choose units in which 7 = 1 . Therefore, the simplest form of the instanton used as a 


trial function in the Sakai-Sugimoto model is 50(4) symmetric. We have generalized this instanton to the deconfined 
geometry in the most natural way, by simply replacing C/kk hr 7 with its analogue in the deconfined geometry U c . 
Since U c is not constant, we cannot work in units where 7 = 1 . It has been argued that corrections beyond the limit 
of infinitely large ’t Hooft coupling A lead to anisotropic instantons which break 50(4) and which might be more 
realistic [45]. For simplicity, we will not allow for these nontrivial configurations which go beyond the BPST ansatz, 
but it is instructive to keep in mind the role of 7 as a parameter for making the instantons anisotropic in a very simple 
way. 

In this paper, we are not interested in a single baryon, but in homogeneous baryonic matter. We thus consider a 
non-interacting gas of Ni instantons [5Q located at positions (X 0n , Z 0n ), n = 1, ...,1V/. We shall assume that all 
instantons sit at the same point in the bulk, at the tip of the flavor branes, Z n 0 = 0 for all n. This is similar to 
the approximation from Ref. [53] . where pointlike baryons were placed at the tip. Our approximation goes beyond 
the pointlike scenario, but does not determine the instanton distribution in the bulk dynamically. Moreover, we shall 
not solve the full X-dependent equations of motion, but rather approximate the instanton distribution by its spatial 
average, such that the locations Xq„ drop out. Therefore, we replace 


4 (p/ 7) 4 if f ,3 x _ 4 (p/7) 4 _ 

K 2 + (p/7) 2 ] 4 v^J ' [(x — x 0 „) 2 + (z/7) 2 + (p/7) 2 ] 4 



JV> 
V ’ 


( 17 ) 
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where V is the 3-volume, and we have defined the normalized function 

3p 4 


q(Z) = 


4(Z 2 + p 2 ) 5 / 2 ’ 


/ oo 

dZq(Z) = 1. 

-oo 


(18) 


Baryon number is generated by the F 2 term that couples to the abelian gauge field Aq in the CS term, because 
the abelian part of the gauge group 17(1) corresponds to the global group at the boundary that is associated with 
baryon number conservation (by choosing the same chemical potential as a boundary value for A 0 at both boundaries, 
Z = +oo and Z = — oo, we ensure that we include a baryon chemical potential, not an axial chemical potential). We 
can thus write the baryon number as |32| 


1 

87T 2 


f d 3 X f dZ Tv[F xj F kz ]e ljk = TV/ , 

J J —OO 


(19) 


where we have inserted Eqs. (15) and ( 0 , re., the baryon number is identical to the number of instantons in our 
instanton gas. We can now write the CS part of the action as 


v n t r°° v f°° 

Scs = -N c — — dZA 0 q(Z) = -JV—ni dza 0 q(z), 


( 20 ) 


with the dimensionless baryon number density ni and the dimensionless function q(z) given in Table [I] [the function 
q(z) depends on the dimensionless baryon width p which, for notational convenience, we shall continue to denote with 
the same symbol as the dimensionful baryon width]. With the help of Eqs. (15) and ( [l7| ), we write all traces over the 
field strengths in terms of the baryon density, 


Tr [i&] 


1 dz niq(u) 


Tr[F, 


(Mkk-R) 3 3 q du ( 2 na ') 2 
where we have introduced 


(M KK i?) 3 

9 u 1 ' 2 


2 y du niq(u) 
3 dz (27ra') 2 

(,Ac ) 2 


dz 

q(u) = 2 q(z') . r-p- / q 3 1 2 ^72 

ou 4v/c ( u 6 — u* + p A u c y > l A 


Tr [FijFku\^ijk 

(MkkR) 3 


/»oo 

/ duq(u) = 1 . 

J u r 


2 niq{u) 

3 ^Trcd ) 2 


( 21 ) 


( 22 ) 


We have included a factor 2 in the definition of q{u) for convenience to ensure tha t q( u) is normalized to one with 
respect to integration over one half of the connected flavor branes. The DBI action (11), with the Tr[E 2 ] terms from 
Eq. (21), together with the CS term (20), yields the action 


V f 00 

s -"tL 


du C , 


with the Lagrangian 


C = C Q 


niq{u) 


iVJt 


IVTcVl + u 3 f T x2 - 
u 2 VJt 


+ 


2 V h 


+ u 3 frx 4 - a; 


- a 0 


where 


£ 0 = + u 3 f T x 4 2 - . 


(23) 


(24) 


(25) 


is the DBI Lagrangian in the absence of instantons. 

Let us discuss the relation of this Lagrangian to the approximation of pointlike instantons. In that case, the 
instanton profile is replaced by a delta function S(u — u c ). The function q(u) is the generalization of that delta 
function. The Lagrangian for pointlike instantons is [53^ 


Tipointlike — T Tlf 


l ^V / /t(u) - a 0 (u) S(u — u c ). 


(26) 


The first term is the same as in our Lagrangian (24), as well as the last term with 5(u — u c ) —> q(u). The second term 
is obtained from the action of TV/ D4-branes that are wrapped on the 4-sphere at u = u c , 


•SYm = NjT 4 j dD 4 dre = N^rn 1 -^- \7/r(n c ), 


(27) 





























FIG. 1: Phase diagram in the plane of temperature and chemical potential in the deconhned geometry with pointlike baryons 
(first computed in Ref. [53]). The baryon onset is a second-order phase transition (thick dashed line), while the chiral phase 
transition is first order (solid line). The thin dashed line would be the chiral phase transition if baryons were ignored (first 
computed in Ref. [63]). At small temperatures, chiral symmetry remains broken for all values of the chemical potential. We 
have rescaled temperature and chemical potential with the asymptotic separation of the flavor branes i. In all plots we use 
dimensionless scales, see Table [T| we do not fit the parameters of the model since in this paper we are only interested in 
qualitative features of the phase diagram, and before the basic features of real-world dense matter are not reproduced by the 
model, any fit would be meaningless. The critical value for tl below which the confined geometry is preferred can be made 
arbitrarily small in the decompactified limit, and thus the deconfinement transition is not seen here. 


with the D4-brane tension X 4 = l/[(27r) 4 f'®]. The energy of the D4-branes i?D 4 = TSd 4 can then be interpreted as 
the mass of Nj baryons. With the help of Table [I] we can write this energy as 

£?d 4 = NjM kk —— - -W/ fr{u c ) • (28) 

47T O 

We recover the baryon mass in the confined geometry with maximally separated flavor branes by setting u c —► «kk = 
4/9 and dropping the temperature-dependent factor. This yields XN c Mkk /(277r) as the mass of a single baryon to 
leading order in A, in accordance with Refs. [32] H2|. In the deconhned geometry, u c \Jfr(uc) /3 depends on chemical 
potential and temperature, giving rise to a medium-dependent baryon mass. We shall come back to this interpretation 
when we compare u c in our calculation with the result from the pointlike scenario, see Fig. [3] 

We show the phase diagram in the plane of temperature and chemical potential for the pointlike approximation in 
Fig. 0 The calculation that leads to this phase diagram was first presented in Ref. |53j . and we recapitulate it in 
appendix [B} Besides the free energy of pointlike baryons, we also compute the free energy of the vacuum (= mesonic 
phase) and the chirally symmetric phase (= quark matter phase) in that appendix. The most important features 
of this phase diagram in our context are the second-order transition from the vacuum to nuclear matter and the 
non-restoration of chiral symmetry for small temperatures and arbitrarily large chemical potential. In the rest of the 
paper we ask the question whether these two unphysical properties can be improved by going beyond the pointlike 
approximation. 


B. All orders in non-abelian field strengths 


We will now discuss a version of the DBI action in which we keep the non-abelian field strengths to all orders. For 
small chemical potentials (which correspond to small baryon densities) the results of this action are identical to the 
ones from the expanded action. For large chemical potentials, the results will obviously differ, and a priori it is not 
clear which version - if any - is “correct”, because there is no unique definition of the non-abelian DBI action to 
all orders in the string tension a'. We shall work with a non-expanded action mainly for two reasons. Firstly, the 
evaluation of the expanded Lagrangian (24) is challenging, and using the non-expanded Lagrangian discussed here 
is considerably simpler, see remarks below Eq. (36 1 . Secondly, we have checked that for the expanded Lagrangian 


there is - within our ansatz - no solution already beyond a moderately large chemical potential, while the solution 
for the non-expanded Lagrangian extends to very large chemical potentials, into a regime where the baryonic phase 
is already disfavored compared to the chirally restored phase. 
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For the non-abelian DBI action we shall use the following prescription. Let us first suppose we work only with 
abelian field strengths Fqu, Fj, Fu ■ Then, we obtain 


det(g + 27 xa'F) = U 4 


R ) 

U 


3/4 


(27ra , ) 2 / T F i 2 £7 + (1 + u 3 f T x 12 - a'o) 


1 + 


(2 W) 2 i* 

2 


(R \ 3 (: 2na') 4 f T (F ij F kU e ijk ) 2 

w 4 


(29) 


Following Ref. [55], we simply replace Fu —> Fj[/ and Fj 
the expressions in Eq. (21), we arrive at the Lagrangian 


Fij and perform the trace of each term separately. Using 


C = u 5/2 s J(1 + gi + u 3 f T x2 - ao 2 )(l + 32 ) - rua 0 {u)q(u ), 


where we have abbreviated 


9i(u) 


f T (u)u 1/2 
3 u 2 c ^ f c (u) 


niq{u ), 


92 (u) 


u 2 c \Jfc(u) 

3u 7 /2 


niq(u). 


(30) 


(31) 


Our goal is to solve the equations of motion for X 4 and do and determine the free energy of the resulting state. 
We work in the grand-canonical ensemble, where the chemical potential is externally given and introduced via the 
boundary value of do, 


a 0 ( 00 ) = ^i. 


(32) 


To be precise, this is the quark chemical potential, such that N c /j. is the baryon chemical potential. The boundary 
condition for X 4 is given by the asymptotic separation L of the flavor branes, 


2 =J d ux ' 4 , 

with the dimensionless separation £ = A'FkL. 

The equations of motion for X 4 and do in integrated form are 

u 5 ^d! oy /l + g 2 


u 5 / 2 u 3 f T x' 4 y/l + g 2 


= n/Q, 


= k , 


•y/l + 31 + u 3 f T x f - a'o 
where k is an integration constant (to be determined), and we have denoted 

u 3/2 yfTc 3 p 2 u c + 2 (u 3 - u 3 c ) 


[ 

Q(u) = / dvq(v) 
J U r 


(u 3 —U 3 + p 2 U c ) 3 / 2 


such that Q(oo) = 1, Q(u c ) = 0. We can solve these equations algebraically for a'o and x' 2 , 

(n/Q) 2 1 + 3i 


an = 


x ' 2 = 


1 + 92 - ^ 57 ^ 
1 + 51 


(»/Q) 2 


u n fr 1 + 32 - 


k 2 

M 8 /t 


1 ( n iQ ) 2 

T „,5 


(33) 


(34a) 

(34b) 

(35) 

(36a) 

(36b) 


At this point, the benefit of using th e no n-expanded Lagrangian (301 becomes apparent: had we worked with the 
expansion in non-abelian gauge fields (24), Oq 2 and x ' 2 would have been solutions to cubic equations, which is much 
more complicated to deal with (unless we had also employed an expansion in the abelian gauge fields). 
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The asymptotic behavior at u = oo of the solutions is 

x '^ u ) = + ’ a 'o( u ) = J/2 +■•• ’ ( 3? ) 

confirming that nj is the baryon density according to the usual AdS/CFT dictionary, which can also be checked later 
numerically by computing the derivative of the free energy with respect to /i. Using these solutions, we can write the 
free energy density U = + S in the useful forms 


n 

Af 



(I + 52 ) 


u 8 f T x ' 4 

k 



- k-ri! 


du 


1 + 32 + 


( niQ ) 


21 


/mi , 


where we have introduced the abbreviation 


(38) 


< = 


' 1 + 91 + u 3 f T x'l - aff u 11 / 2 f T x ' 4 u 5 / 2 d ' 0 


V 1 +9i 


1+32 


njQ 




1 +32 - 


fc 2 

u s f T 


( niQ ) 2 
7/5 


(39) 


This free energy is divergent at the holographic boundary u —► 00 . Replacing the upper boundary of the integral by 
a cutoff A, one finds that the divergent contribution is 2/7 A 7 / 2 . All phases we consider in this paper have the same 
divergence, which is a pure vacuum contribution, i.e., does not depend on /i or T. Therefore, and since we are only 
interested in differences of free energies, we can simply drop this contribution. 

Besides the functions x 4 (u ) and &o(u), the system contains the parameters nj, p 1 and u c , and we have to minimize 
the free energy with respect to them. At first sight it seems curious to minimize with respect to the density n/, since 
we work in the grand-canonical ensemble where the density should be given as a function of p. However, in our setup 
nj should primarily be considered as a parameter of the Lagrangian. This parameter turns out to be identical to the 
baryon density. This is not necessarily the case, since the total density may receive contributions from other sources, 
for instance through a magnetic field [ 53j . 

In order to minimize the free energy with respect to nj, p , and u c , it is convenient to read the Lagrangian as a 
functional of ao(u),a' 0 (u),x' 4 (u),ni, p , u c , where the three functions do(u),a' 0 (u),x' 4 (u) also depend on nj, p, u c . Then, 
the minimization with respect to nj can be written as 


(dc^d do r 00 au 

V da' 0 dn / + dx' A dm ) M=Uc + J Uc U dm ’ 


(40) 


where we have used the equations of motion. All boundary terms vanish because the derivative is taken at fixed 
p = fio(oo) and fixed separation 1/2 = £ 4 ( 00 ) — 2 : 4 (u c ), and because 


dC 
da 0 



(41) 


due to Eq. (34a). Consequently, only the explicit derivative with respect to nj remains. The derivative with respect 
to the baryon width p is taken completely analogously. The derivative with respect to u c is a bit more subtle. We 
obtain 


(dt^dd 0 dC^dx^y- 00 _ , r°° dCy 

\da ' 0 du c + dx ' 4 du c ) u=Uc v ~ v - + J Uc U du c ' 


(42) 


Now we need to take into account that do(u c ) and x 4 (u c ) depend on u c explicitly and through the dependence on u. 
In the boundary terms, however, only the explicit dependence is relevant. We thus write 


dx 4 
du c 



dx 4 (u c ) 

du c 


x' 4 (u c ) 


(43) 


(and the same for So), where the left-hand side is needed in Eq. (42), the first term on the right-hand side denotes 
the full dependence on u c , and the second term the dependence via u. Using this relation and 


dC 

dx ' 4 


= k , 


(44) 
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which follows from Eq. (34b), it turns out that the boundary term at u = u c gives a contribution kx' 4 (u c ), i.e., Eq. 
(42) becomes 


r c 

0 = (kx 4 - C) u =u c + / 

J U r 


du 


dC 

8Un ' 


In summary, the three equations to minimize the free energy can be written as 


0 = 


0 = 


du 


du 


«' 5/2 (dgif-i , ^ 


L 2 V5n / C ” 1 + £ C)+ “°° 


5/2 (dg,_ i , dg 2 


C ) + n I^O~a~ 


dp 


dp 


.,dQ 


dp 


(45) 

(46a) 

(46b) 


0 = -u 5 / 2 J [1 + gi (u c )\ 


l + g 2 {u c ) - 


k 2 

ulfr{u c )_ 


du 


v 5/2 (d9i c -1 dte 

du c du c 


, SQ 


(46c) 


In all three equations, we have eliminated do in favor of dp via partial integration. This is advantageous because 
do can only be obtained by numerically integrating Eq. (36a), i.e., by eliminating do we avoid two nested numerical 
integrations. 

We have arrived at a system of coupled algebraic equations - Eqs. (46) plus the equation for the asymptotic 
separation ( |33| ) - which has to be solved for k, m, u c , and p for given chemical potential p and temperature T [which 
appears in fr(u)}. Before discussing the results, let us comment on the numerical evaluation of these equations. 


C. Numerical evaluation 


The equation that requires some explanation is the minimization with respect to u c (46c), for which we need to 


know the behavior of various functions at the tip of the connected flavor branes, u — u c . We find 


niq(u) = 


yjU U f 


+ 0[(u — Uc ) 1 / 2 ], niQ(u ) = 2\/3 au 2 c \/u — u c + 0[(u — u c ) 3 ^ 2 ], 


(47) 


which implies 


9i(u) 


af T (u c )u c 
3(u - u c ) 


+ 0 ( 1 ), 


g 2 (u) = a + 0(u — u c ), 


(48) 


where we have abbreviated 


3 nj 


Apu, 


3/2 ' 


Inserting these expansions into the solution for x' A ( 36b), we find that x' 4 diverges at u = u c 

Ci . , nr / , 1 / 9 , 1 a 1 / 2 k 


x\(u) = 


y/u - li ( 


+ 0[(u - u c ) 1/2 }, Cl = 


VSu c y/ulf T (u c )(l + a) - k 2 


(49) 


(50) 


Consequently, the embedding of the flavor branes x 4 (u) is smooth at u = u c for all values of nj. This shows that 
the cusp, introduced by the approxi matio n of the delta-like baryons, is, not surprisingly, removed by the smooth 
instantonic baryons; see also appendix B 1 where we review the pointlike approximation and see that in that case x' 4 
is finite at u = u c and only diverges for n/ —> 0. 

We can now compute 


1 

k 2 

[1 +gi(u c )] 

l + g 2 (u c ) 8 , , 

U°fT(Uc) J 


ak 


3u 2 Ci^/u — u c 


+ 0[(u - U c ) 1/2 ] 


(51) 
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x ' 4 

q,p, ci 

X 4 

Q, 00 , 51 , 32 , a, C 

u,u T , p,a 0 , p 

ni 

n/N 

k 

e 


i 

1 

r 2 

£~ h 

r 7 

£ -s 

-3/2 

U c 

U c 

- 1/2 

U c 

1 

U c 

u 5/2 

Uic 

7/2 

Uc 

ui 


TABLE II: The system of four coupled equations ( |33| ), ( |46a| ) - ( |46c| is most conveniently solved by first rescaling all quantities 
with the constant £, and then with the variable u c , as given in this table. As a consequence, £ disappears from the equations, 
and u c only remains in Eq. (331, which thus decouples from the other three equations. 


i.e., we have obtained a divergent contribution (and no constant term). However, the integral in Eq. (46c) also contains 
a divergent term, and both divergences exactly cancel each other. The divergent term of the integral arises in 


dgi x dg 2 
du c du c 


9i(u) 


P(u ) 


2 

U r . 


C iu)+g 2 {u) 


p(u) 



gi{u) df c 
fc(u) du c 


C *(“)> 


(52) 


where we have abbreviated the function 


2 5 3 u 2 — p 2 

p(u) =- 1 - — —3 - 3 2 — ’ 

u c Zw 5 — u% + p u c 

which originates from taking the derivative of q(u) with respect to u c . Now, with 

dQ u 3 + 2u 3 


(53) 


(54) 


we can write Eq. (46c) as 


0 = 


du 


5/2 


C<? 2 (p+^) +C 


u 3 + 2 u 3 

niqa 0 


ak 


3 u 2 c C Si 


3 u 2 u c §u\c\ (u — u c ) 3 / 2 u x / 2 f c 


(55) 


and one easily checks that the divergences in the two last terms cancel each other, rendering the integral finite. The 
other equations (33), (46a), and (46b I, are evaluated much more straightforwardly and require no further explanation. 


For the numerical evaluation, we first note that rescaling of all quantities with appropriate powers of i eliminates i 
from all equations. The various powers of i are given in Table [TT| We shall not introduce new symbols for the quantities 
rescaled with £, but rather write l explicitly in all plots. Then, we can proceed analogously with u c (remember that l 
is a constant parameter, while u c is a dynamically determined variable). Firstly, this removes the variable u c from the 
lower boundary of the integrals, which is convenient for the numerical evaluation. And secondly, after this rescaling, 
u c drops completely out of Eqs. (46a) - (46c), i.e., we can solve them for (the rescaled) k, ni, p and afterwards 
compute u c from Eq. (33), which is then used to undo the rescaling of k. nj , p. 


We also note that p only appears in Eq. (46a), and thus one can further reduce the number of coupled equations to 
two if one fixes nj and solves for k and p 1 and afterwards computes the resulting p and u c . Finally, we mention that 


the numerical integration is best done after a change of variables according to Eq. (131 from u to z wit h z defined as 
u = (1 + z 2 ) 1 / 3 (after scaling out u c ). This is particularly helpful for the numerica 
the most challenging one. 


integration in Eq. (55), which is 


D. Second-order baryon onset and chiral restoration 


In this section, we present our numerical results for the instanton gas. Our emphasis is on the deconfined geometry. 
Nevertheless, we start with the confined geometry with maximally separated flavor branes. This allows us to compare 
our results to the literature and serves as a warm-up exercise for the more difficult calculations in the deconfined 
geometry. All relevant equations for the confined geometry are collected in appendix ED the calculation basically 
reduces to solving Eqs. (A20), the stationarity equations for the baryon density n/ and the baryon width p. In contrast 
to the deconfined geometry, the results do not depend on temperature. The location of the tip of the connected flavor 
branes is fixed in the case of maximally separated branes and identical to ukk- Baryon width, baryon density, Oo(ukk), 
and the free energy density are plotted as a function of the chemical potential in Fig. [2] We see that for small n/ we 
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Confined geometry, instanton gas 





FIG. 2: Confined geometry with maximally separated flavor branes: baryon width p (upper left panel), baryon density ni (upper 
right panel), value of ao at the tip of the connected flavor branes (lower left panel), and free energy density (lower right panel), 
all as a function of the chemical potential p. (Note the larger p scale in the lower right panel.) The dashed lines correspond to 
the results for pointlike baryons, where p = 0. The plots show the second-order baryon onset at p = ukk/3 = 4/27 ~ 0.148. 
Below the onset, i.e., in the mesonic phase, oo(mkk) = p, ni = 0, Q. = 0, and p undefined. For the dimensionless units used 
here and in all following plots, see Table [T] 


reproduce the results of the pointlike approximation, in particular we find a second-order baryon onset 3 . The onset 
occurs at p = wkk/ 3 = 4/27. With the help of Table [I] and the comments below Eq. (281 we see that this critical 
chemical potential is identical to the baryon mass, which implies that there is no binding energy. The plot of the 
free energy shows that at very large chemical potentials the finite-size baryons become energetically more costly than 
pointlike baryons. At even higher chemical potentials, p > 2.37, beyond the scale shown in Fig. [2j we do not find any 
solution. (There is a second solution up to that chemical potential p ~ 2.37 which we have not included in the plots 
since its free energy is always higher than the one of the solution plotted here). 

We now turn to the deconfined geometry, whose results are shown in Figs. [3] (zero temperature) and [4] (phase 
diagram). In the upper left panel of Fig. [3] we show the baryon density ni and compare it to that of the pointlike 


3 In Ref. I56| it was claimed in an almost identical calculation that the baryon onset becomes first order within the approximation of the 
instanton gas. Even though in that reference the DBI action was expanded for small non-abelian gauge fields, the transition to baryonic 
matter should be of the same order because our result reduces to that of the expanded action for small baryon densities. However, in 
Ref. [56| . ao(^KK) has been set to zero. In doing so, a lower boundary for the baryon density is “forced” upon the system (and the 
free energy is minimized only with respect to p, not also with respect to nj). Our calculation shows that if uo(^kk) is determined 
dynamically, it chooses to approach p at the baryon onset, connecting continuously to the mesonic phase, where do(u) = p for all u , see 
lower left panel of Fig. [2] As a consequence, the baryon density vanishes at that point and the onset is second order. 
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Deconfined geometry, instanton gas, T = 0 






0 20 40 60 80 100 


jjl ' 2 


FIG. 3: Deconfined geometry at zero temperature: baryon density ni (upper left panel), baryon width p (upper right panel), 
location of the tip of the connected flavor branes u c (lower left panel), and ratio of the pressures of baryonic and quark matter 
phases (lower right panel), all as a function of p and rescaled with the asymptotic separation of the flavor branes l. (Note 
the larger p scale in the two lower panels logarithmic in the lower left panel.) The dashed lines are the results for pointlike 
baryons. Due to numerical difficulties for very small p, the solid lines do not connect with the lines of the mesonic phase below 
the second-order baryon onset. In the lower right panel, the solid line stops at a point beyond which we have not found any 
solutions (and where baryonic matter is already disfavored). 


approximation. As for the confined geometry, we observe a second-order phase transition whose properties around 
the transition (i.e., for small ni ) are well approximated by pointlike baryons. This transition is reminiscent of a 
Bose-Einstein condensation [55] . which is maybe not surprising since at infinitely large N c baryons may lose their 
N c = 3 fermionic nature. Our numerical result for the instanton gas is not shown all the way down to the onset 
because the numerical evaluation becomes problematic at very small baryon densities. This regime corresponds to 
very small baryon widths p , as the upper right panel shows, and thus we have to integrate numerically over highly 
peaked functions, which is hard. However, the plots show that our calculation works well down to a regime where the 
result already merges with the pointlike approximation. And, in the confined geometry, where we can perform the 
calculation down to the onset, we have seen that instanton gas and pointlike approximation become indeed identical. 
Therefore, it would be very surprising if the full result in the deconfined case did not follow this approximation down 
to the baryon onset, even though our numerics cannot prove it. 

The lower left panel shows the location of the tip of the connected flavor branes u c . This is an interesting quantity 
because it can be related to the medium-dependent baryon mass, at least in the pointlike approximation, see Eq. (28) 
and comments below that equation. Interestingly, the behavior of u c for the instanton gas at large chemical potentials 
differs qualitatively from that of the pointlike approximation. While instantonic baryons appear to become heavier 
without bound, the pointlike baryon mass saturates at a finite value. As a consequence, one can expect the instanton 
gas to become energetically more costly with increasing p relative to the gas of pointlike baryons. This expectation is 
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FIG. 4: Phase diagram in the deconfined geometry with baryons in the instanton gas approximation. Since for small baryon 
density the result approaches the one from pointlike baryons, the second-order baryon onset (as well as the transition from 
the mesonic to the chirally restored phase) are identical to the ones in Fig. [I] The first-order transition between baryonic and 
chirally symmetric phases is shown by the thin (blue) line for pointlike baryons and by the thick line for instantonic baryons. 
The main qualitative difference is the chiral restoration for small temperatures at very large fx (note the logarithmic /x scale, 
while the /x scale in Fig. [T| is linear). The black line for the instanton gas is not connected to the critical point where the 
second-order onset line meets the first-order chiral transition line because of numerical difficulties at small baryon width, see 
also Fig. [3] 


borne out in the lower right panel: the pressure P = —0 of pointlike baryons approaches that of the chirally restored 
phase, P oc /x 7 / 2 , from above thus being favored over quark matter for arbitrarily large /x. (The free energy 
of the chirally restored phase is derived in appendix B3 ) In contrast, while the instanton gas is energetically less 


costly (= has higher pressure) at small /x than the pointlike baryons, the curve of its pressure crosses the pressure 
curves of the pointlike approximation and the one of chirally restored quark matter. This shows that in the instanton 
gas approximation there is a zero-temperature chiral phase transition to quark matter at large chemical potential. 
The corresponding critical chemical potential is extremely large, about 200 times larger than the chemical potential 
for the baryon onset (/x ~ 35/T 2 , while the onset occurs at /z ~ 0.17/T 2 ). Therefore, if we were to fit the baryon 
onset to QCD (which makes little sense due to its obviously unrealistic properties), our calculation would predict 
chiral restoration to occur at a baryon chemical potential of about 200 GeV. We have to keep in mind, however, 
that the result at large /x has to be taken with a lot of care because of various reasons: we have used a simple, not 
uniquely defined, prescription of the non-abelian DBI action, whose validity is guaranteed only for sufficiently small 
/x; we have employed the probe brane approximation, which becomes questionable for very large gauge fields on the 
flavor branes; and, the instantons start to overlap at densities well below the critical chemical potential for the chiral 
phase transition, which challenges our instanton gas approximation. For these reasons we do not attempt to use this 
approximation for any quantitative predictions, since clearly it first has to be improved in the future. Nevertheless, 
the existence of chiral restoration of baryonic matter in the Sakai-Sugimoto model at high densities, at least within a 
certain approximation, is a very interesting and - to our knowledge new observation and important progress towards 
a high-density, strongly coupled model outlined in the introduction. 

We can now compute the phase transition between baryonic and chirally restored phases for all temperatures. The 
resulting phase diagram is shown in Fig. [4] 


IV. BARYONIC MATTER FROM A HOMOGENEOUS ANSATZ 
A. Ansatz and calculation 


So far, we have worked with an ansatz for the non-abelian gauge fields that was based on the BPST instanton 
solution. In this section, we discuss a second ansatz, given by (5S| 


A( z ) = -<*% 


H{Z) 


A-z = 0 , 


2 


(56) 
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with some function H(Z) 1 to be determined. With + i[A^, A„\, the field strengths become 


Fj j j /, (71,' 


H 2 


F iZ = <Ji 


dzH 


(57) 


where [a,. a 7 \ = 2 ieijkcrk has been used. This ansatz is different from the instanton solution (A 8 ), where both Az 
and Ai are nonvanishing and depend not only on the holographic coordinate Z but also on position space. In our 
approximation of the previous section, we have integrated out the position space dependence before solving the 
equations of motion, i.e., eventually also our instanton gas was spatially homogeneous. Nevertheless, we shall refer 


to Eq. (56) as “homogeneous” since it does not involve any spatial dependence to begin with. As noticed already in 
Ref. |58|. the ansatz (56) only yields a nonzero baryon density if A^Z) is allowed to become discontinuous at Z = 0. 
The reason is that the baryon density assumes the form 


1 r°° If 00 


(58) 


and H(±oo) = 0 has to be required in order to ensure a finite free energy (up to a constant divergent term that 
we always subtract). Therefore, if H 3 was a continuous function, the integral (58) would simply give H(+ oo ) 3 — 
H{— oo ) 3 = 0. We shall briefly return to this issue below Eq. (70), when we know the exact behavior of H(Z) around 
Z = 0. 

As for the instanton gas, we shall work in the coordinate U and on one half of the connected flavor branes [for the 
relation between Z and U see Eq. With Tr [crjcq] = 6 we have 


Tr [Ffu] = 


3h' 2 


Tr[F 2 


3A 2 0 h 4 


Tr [FijF kU ]eijk _ 3A 0 h 2 h! 


2(27ra0 2 ’ ( M KK R) e {2na') 2 ’ ( M KK R ) 3 (27ra') 2 ’ 

where the dimensionless gauge field h is defined analogously to do, see Table |T[ and where we have abbreviated 

\ _ A 
Ao= 4^' 


(59) 


(60) 


Replacing all abelian F 2 terms in Eq. (291 with these non-abelian versions and using the CS term (| 8 |) yields the action 

3A 2 d 4 A _ „ . 9 h 2 h r 


/»oo 

[7 

S = J\f / du 

J u c 

" 5 , V( 


3 f T h’ 2 3 /2 „ 

1 H---h u 3 f T x'i - a j 


1 + 


2 u 3 J 


clqXo- 


(61) 


In contrast to the instanton gas, the Lagrangian now depends on the’t Hooft coupling explicitly, there is no rescaling 
of the fields by which we can get rid of A. The reason is the asymmetric ansatz for Az and A,;, which leads to a 
different scaling behaviour of Fij and Fuj. Even though this ansatz is somewhat more simplistic than the instanton 
gas, the explicit appearance of A allows us to capture some physics away from the A = oo limit. We shall see that 
only at finite A we obtain a realistic first-order onset of baryonic matter. Of course, since we do not go beyond the 
classical gravity approximation, we cannot claim to include finite-A corrections systematically and have to consider 
our results at finite A as an extrapolation. 

The Lagrangian now contains the additional unknown function h(u), besides do(u) and 2 : 4 ( 14 ). We thus have to 
solve an additional Euler-Lagrange equation. This seems to render the calculation more tedious, but we shall see 
that the three differential equations can be decoupled, and eventually the effort needed to evaluate this approach 
numerically is comparable to the one needed for the instanton gas. Employing the same compact notation as in Sec. 
|IIIB| we abbreviate 


3 f T h'(u) 2 
9i{u) = -x- 


_ 3A lh{uY 

*<“> = 


The baryon density is obtained from the CS term, 

3A 0 f 

n i = „ / 


> O \ 

dud u h 3 = — ^h(u c ) 3 , 


(62) 


(63) 


where we have used h{ 00 ) = 0. (For the sake of a consistent notation, we continue to denote the baryon density by 
n/, although the current approach is not based on the instanton solution.) 
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The equations of motion for do and X 4 are 


J - 5// ~a'oVl + 92 


V 1+ 9i+ u 3 f T x'i - a' 0 

U 5/2 U 3 f T X 4-y/l + ff 2 

\/l + 9i + u 3 f T x-2 - d'o 


= n/Q. 


= fc, 


with the integration constant k and, in analogy to the function Q from Eq. (35) 


Q{u) = f^[M u ) 3 - h(u c ) 3 ] = 1 + ^-h(u ) 3 . 

2,ni Zni 


(64a) 


(64b) 


(65) 


(To emphasize the analogy to the calculation of Sec. Ill we use the same symbols g \, g 2 , Q , but of course need to 
keep in mind that they denote different functions.) In complete analogy to Eqs. (36) we thus find 
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x'i = 


l + 3i 


ullf Tl + 92-^ + ^' 


For the equation of motion for h we need 
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dh 
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where 


c = 
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(67) 


( 68 ) 


The expression on the right-hand side does not contain do and X 4 anymore, which can thus be eliminated from the 
equation of motion for d, 




_idgi\ 9\ 0 h 2 (mQ _ 5/2 dg 2 
u-5/2 “ u ** dh ' 


(69) 


We are thus left with this single differential equation for h which has to be solved numerically. From Eq. (63) we 
know the boundary value at u = u c in terms of nj. By expanding Eq. (691 around u = u c , we find the behavior 


M«> = - Ud 


2»,y / y .— 

— + aiy/u - u c 

0 A 0 J 


+ a 2 (u - u c ) + ... , 


(70) 


with coefficients ai, a 2 . Inserting this expansion into Eq. (66b), we see that x ' 4 oc 1 /y/u — u c for u close to u c . This 
divergent derivative shows that the embedding of the flavor branes is smooth, just like for the instanton gas, i.e., 
there is no cusp at the tip of the connected flavor branes, which occurs in the pointlike approximation. Remember 
from the discussion at the beginning of this section that the function h(z) must be discontinuous at the tip. Thus, 
if we want to continue h(z) into the second half of the connected flavor branes, we have to add a factor sgn(z), such 
that h is antisymmetric under z —► —z. As a consequence, d z h 3 is symmetric and continuous, and the z < 0 half of 
the connected branes gives the same contribution to the baryon density as the z > 0 half. The function d z h 3 is the 
analogue of the instanton profile q(z), see Eq. (18). In contrast to the instanton profile, d z h 3 has a cusp at z = 0. But 
unlike the delta-like baryons from the pointlike approximation, d z h 3 is finite at z = 0 and has a finite width. Having 
this picture of the complete connected flavor branes in mind, we can now go back to working on one single half. 


For a given a±, we can compute k: from the requirement that Eq. (69) be fulfilled to lowest order, (u — u c )°, we 


find 


k 2 = ulf T {u c )- 
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8 u c + 9 a\f T {u c ) 


(71) 
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(The higher-order terms contain the higher-order coefficients 02 etc., which can be computed successively as a function 
of a\ in this way, but will not be needed in the following.) We can thus eliminate k in favor of a\. For a given baryon 
density nj, which determines the boundary value h(u c ), we determine the coefficient a\ numerically via the shooting 
method, such that h(o o) = 0 . 

Next, we need to take into account the minimization of the free energy with respect to n/. Reading the Lagrangian 
as a functional of do, d' 0 , x' 4 , h, h! (with no further explicit dependence on ni), we obtain 


0 = 


dC ddo 
da' 0 dni 


U — OO r\ /1 r\ 

ujL ox 4 

u -°° dC dh 

U—OO 

—K 

dC 

u=u c + dx 'i dn I 

U=u a + dh ' dn I 

- 7 

U—U c 

dh' 


= 0, 


(72) 


where we have used that &o(oo) = p, £ 4 ( 00 ) — £ 4 (u c ) = 1 /2, h( 00 ) = 0 are fixed and thus do not depend on n/, and 
I — 0 — I — k 

da' I“=“c — dx' I»=»= — 


With the help of Eq. ( |67| , the condition ( |72j ) can be used to compute do (u c ), 

2 u 8 


do [u c ) = 


u c \J fr{u c ) 


3A lh 4 (u c ) 


1 - 


k 2 


Ifriuc) 


(73) 


From this relation we can already see that for large ’t Hooft couplings, Ao —> 00 , we recover the result of pointlike 
baryons, see Eq. (B 6 ). We shall come back to this observation in our discussion of the numerical results. 


Recall that in the case of the instanton gas we additionally had to minimize with respect to the instanton width 
p and the location of the tip of the connected flavor branes u c . Now, there is no parameter p , and the minimization 
with respect to u c is automatically included in the equations of motion, which can be seen as follows. Minimizing the 
free energy with respect to u c yields 
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dC ddo 
da' 0 du c 
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dC dx 4 
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where we used the same arguments as explained below Eq. (|42|). We compute 
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Since Q{u c ) = 0, we obtain 


kx ' A + h 'dh' c 


u 5 / 2 k 2 

V 1 + gi V + 92 u 8 f T 


= 0, 


(74) 


(75) 


(76) 


because gi oc (u — u c ) 1 , and 52 = const at u = u c , i.e., we are already sitting on a stationary point of the free energy 
with respect to u c once we fulfill the equations of motion. 


Finally, we need to fulfill the boundary conditions (32) and (33). For a given chemical potential p, these two bound¬ 
ary conditions, together with the equation for h (69), determine a\, ni, and u c . Instead of solving this complicated 


system of equations simultaneously, we can proceed as follows. First, we rescale all quantities by appropriate powers 
of it c and the asymptotic separation £, as explained in Sec. |III C[ see Table [ll] In addition to the quantities listed 
in that table, we need A ~ u c 1 ^ 2 , ai ~ v/J 2 , h ~ u c and accordingly for using i ~ u~ 2 . Then, t is completely 
eliminated from the equations, while u c 


only appears in a trivial way in the boundary condition (331, which thus 


decouples. We now choose a value for (the rescaled) n/, and determine ui and h{u) via the shooting method from 
Eq. (69). Then, after determining do{u c ) from Eq. (73), we use the boundary condition (32) in the form 


9 = 


dn ( 


do{u c ) 


(77) 


to compute (the rescaled) p. [Computing p as a function of n/ is advantageous also because it is a single-valued 
function, while n/(/x) is two-valued.] The disadvantage of this procedure is that we only work with rescaled quantities, 
and the rescaling with u c has to be undone to obtain the final results (the rescaling with t is trivial because £ is a 
constant). As a consequence, we have used this procedure for the zero-temperature phase diagram in Fig. [7J but it 
cannot be used when we wish to work with a fixed’t Hooft coupling A. In this case, we must not rescale A, and the 
numerics become somewhat more difficult because we now have to simultaneously solve the differential equation for 
h and the boundary condition (33) to determine u c . This has been done for the plots in Fig. [ 6 ] 


For the free energy comparison to the mesonic and chirally symmetric phases, we notice that, due to our choice of 
notation, the free energy has exactly the same form as given in Eq. (38), only with different functions £, 52 , and Q. 
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Confined geometry, homogeneous ansatz 



A 


A 


FIG. 5: Upper panels: first-order baryon onset for the confined geometry with maximally separated flavor branes within the 
homogeneous ansatz at a fixed ’t Hooft coupling A = 10. The left panel shows the baryon density, which is discontinuous at 
the onset, while the right panel shows the free energy which favors baryonic matter for S7 < 0. In both panels, the dashed 
continuation of the curve shows the metastable solution (up to the turning point) and the unstable solution. Lower panels: 
baryon density just above and chemical potential at the first-order baryon onset as a function of the’t Hooft coupling A; the 
critical chemical potential approaches that of the second-order onset in the pointlike approximation fi = ukk/3 for A —> oo. 


B. First-order baryon onset and \-fi phase diagram 


Again we start our discussion of the numerical results with the confined geometry with maximally separated flavor 
branes, whose equations for the homogeneous ansatz are collected in appendix |A 2| The results are shown in Fig. [5] 
The two upper panels, where we plot the baryon density and the free energy as a function of /z at a fixed value of 
A, show a first-order baryon onset. This is in accordance with Ref. [58] . where this observation was already made 
for the confined geometry, but without a full numerical evaluation. The lower two panels show the dependence on 
the’t Hooft coupling A of the baryon density just above the onset ro/ i0nse t and of the chemical potential at the onset 
Monset - We see that n/ i0n set vanishes for A —> oo, while /z onset approaches that of the second-order onset of the pointlike 


approximation and the instanton gas of Sec. Ill 


In the deconfined geometry we restrict ourselves to zero temperature and present our results in Figs. [6] and [7] In 
the upper left panel of Fig. [G] we plot the baryon density as a function of the chemical potential for three values of the 
’t Hooft coupling, \/(. = 20,40,10 3 . As for the confined geometry, we obtain two solutions, but only show the stable 
branch to keep the plot simple. We see that the jump in nj becomes small for large A (although the size of the jump is 
not a monotonic function of A), and the critical chemical potential approaches that of the pointlike approximation. In 
the upper right panel we show the free energy of the solutions, in comparison to the mesonic and the chirally restored 
phases. Here we have included the metastable and unstable branches, whose free energy is larger than that of the 
mesonic phase for all m- In the lower left panel we show the location of the tip of the connected D8- and D8-branes 
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Deconfined geometry, homogeneous ansatz, T = 0 



FIG. 6: Upper left panel: zero-temperature baryon density for three different values of the ’t Hooft coupling as a function of 
the chemical potential, showing the discontinuity at the first-order transition. In all panels of this figure, the dashed line is 
the result for pointlike baryons. Upper right panel: corresponding free energy, including the metastable and unstable branches. 
The thin dashed-dotted lines are the free energies of the mesonic phase (horizontal line) and the chirally restored phase. Lower 
left panel: location of the tip of the connected flavor branes on a large, logarithmic p scale. Lower right panel: ratio of the 
pressures of the baryonic and chirally restored phases for two large values of the’t Hooft coupling on a large, linear p scale. 


on a very large p scale. We see that u c saturates at a finite value for p —► oo, just like for pointlike baryons and in 
contrast to the instanton gas approximation. The value at infinitely large A and p depends on the order in which these 
limits are taken: if p is fixed to any arbitrarily large, but finite value, u c approaches the pointlike result, u c — 1.69/£ 2 , 
for A —> oo; if, on the other hand, A is held fixed, u c approaches a smaller value for p —» oo, namely u c — 1.23/£ 2 . 
Finally, in the lower right panel, we show the ratio of the pressures of the baryonic and chirally restored phases, in 
analogy to the lower right panel of Fig. [3j This plot demonstrates that there is no chiral restoration: for any finite 
value of A, the ratio approaches a finite value larger than 1 for p —> oo. 

We may now compute the critical chemical potential for the baryon onset for all A. The resulting zero-temperature 
phase diagram in the A -p plane is shown in Fig. [7] As the results of Fig.[6]have already suggested, we recover the result 
of the pointlike approximation at very large A: the baryon onset becomes a weak first-order transition (second order 
for A -» oo), and we even reproduce asymptotically - within the numerical accuracy and performing the calculation 
up to A ~ 1.4 x 10 4 - the numerical value for the critical chemical potential, p ~ 0.17/£ 2 . Going to smaller values of A, 
the critical chemical potential increases and the system behaves as just discussed and shown in the other three panels. 
For sufficiently low values of A, the situation differs qualitatively. In this case, the mesonic phase is not superseded 
by baryonic matter, but a chiral phase transition to quark matter occurs, and only at very large p baryonic matter 
becomes favored, resulting in another chiral phase transition. The topology of the phase structure in the A -p plane 
gives rise to a tricritical point at (A ,p) ~ (10.1£, 0.44/£ 2 ), where all three phases coexist. 

The main conclusion from this phase diagram is that, at fixed ’t Hooft coupling, once baryonic matter is created in 
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FIG. 7: Phases in the plane of chemical potential g and’t Hooft coupling A at zero temperature for the homogeneous ansatz 
of the non-abelian gauge fields. All solid lines are first-order phase transitions. The thin dashed line marks the (second-order, 
A-independent) baryon onset for the instanton gas. For not too small values of the ’t Hooft coupling, there is a first-order 
phase transition from the vacuum to baryonic matter, but - if the coupling is kept fixed - no subsequent transition to chirally 
restored matter. 


a realistic first-order phase transition, it never wants to disappear again, i.e., the baryon onset is never superseded by 
chiral restoration at high densities, as one would expect from QCD. In this sense, the homogeneous ansatz behaves 
exactly opposite to the instanton gas, where the onset was an unrealistic second-order phase transition, but there was 
chiral restoration at high densities. 

Interestingly, the phase diagram of Fig. [7] shows that if the ’t Hooft coupling were allowed to vary as a function 
of /lz, it would be possible to move from the vacuum through a first-order onset into the baryonic phase and then 
via a chiral phase transition to the quark matter phase, as expected from QCD. To this end, the ’t Hooft coupling 
A = g 2 N c would have to decrease with g, in accordance with the running of the QCD coupling g. It is thus tempting 
to speculate what the trajectory of the running coupling would be in the phase diagram of Fig. [7] Suppose we reach 
the quark matter phase with a running coupling. Then, we know that at some sufficiently large g the QCD coupling 
runs logarithmically, g 2 oc 1/ln/i. The transition line between the quark matter and baryonic phases in Fig. ul 
however, decreases with A through a power law (the numerics suggest g oc A -5 / 2 , but computing this phase transition 
line numerically becomes increasingly difficult for large g). As a consequence, it seems the system necessarily has to 
re-enter the baryonic phase at large g. But, we need to keep in mind that the current model cannot be expected 
to be valid in the asymptotically dense regime where the system becomes weakly coupled. Therefore, it is perfectly 
conceivable that the (strong) coupling runs through the three phases of Fig. [7] in the “right order” before we anyway 
have to stop trusting our calculation. The chiral phase transition is expected to occur at moderate, not asymptotically 
large, densities. Of course, in the present approximation, there is no prediction for how A might run, and thus here we 
do not make any attempts to model such a running. Possibly this question can be addressed by taking into account 
the backreaction of the flavor branes on the background geometry, which appears to result in a running coupling, 
although different from QCD [38]. 


V. CONCLUSIONS 

We have studied cold and dense baryonic matter in the Sakai-Sugimoto model and have addressed the question 
whether the model can account for a first-order baryon onset and a zero-temperature chiral phase transition from 
baryonic to quark matter. To this end, we have focused on the deconfined geometry and the decompactified limit, where 
the separation of the flavor branes at the holographic boundary is small compared to the radius of the compactified 
extra dimension of the model. To our knowledge, this is the first time that holographic baryonic matter beyond the 
simple pointlike approximation has been considered in this setup. As a reference calculation and for a comparison 
to existing results in the literature we have also considered the confined geometry with maximally separated flavor 
branes. 

In our description of baryonic matter we followed two different approaches. Firstly, we have considered an instanton 
gas [5B], which makes use of the BPST instanton solution that has been employed to study vacuum properties of 
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baryons in the Sakai-Sugimoto model (42] . This approach can be considered as a generalization of the pointlike 
approximation for baryonic matter !53! . We have recovered the pointlike approximation in our calculation for the case 
of small baryon densities. It has turned out that the baryon density is allowed to become arbitrarily small, resulting 
in a second-order phase transition from the vacuum to baryonic matter. At large densities, however, the instanton 
gas differs qualitatively from pointlike baryons. Most notably, the medium-dependent baryon mass increases without 
bound, such that baryonic matter becomes energetically more and more costly. As a consequence, and in contrast to 
the pointlike approximation, chiral symmetry is restored at zero temperature at (extremely) large densities. 

Secondly, we have employed a homogeneous ansatz for the non-abelian gauge fields, which only depends on the 
holographic direction |58| . This approach, in contrast to the instanton gas approximation, allows for a nontrivial 
dependence on the ’t Hooft coupling A and thus for going away from the infinite coupling limit (at least in an 
extrapolating sense since we do not compute finite-A corrections systematically). This is important because only for 
finite values of the’t Hooft coupling do we see a first-order baryon onset, as expected from QCD. If A is kept fixed and 
not too small, we have found that after baryonic matter is created in a first-order phase transition, it remains favored 
for all densities, i.e., there is no chiral restoration. A sequence of phases with increasing chemical potential that is 
expected from QCD (vacuum - nuclear matter - quark matter) is only possible if A were to decrease with increasing 
chemical potential. Although this is an intriguing observation since this is exactly what is expected from the QCD 
coupling, the present approximation does not include any running of A. The main conclusion from our two different 
approaches is thus that neither one shows a first-order baryon onset and chiral restoration at high density, but we see 
chiral restoration in the instanton gas approximation and a first-order onset in the homogeneous approach. 

In order to interpret this conclusion we need to keep in mind that we have employed several approximations and 
extrapolations. We have used a simple prescription for the non-abelian DBI action and neglected any backreaction of 
the flavor branes on the metric. Therefore, our results for very large densities have to be taken with some care. We 
have used the simple flat-space, ^0(4) symmetric BPST instanton solution. In going from a single instanton profile 
to a gas of instantons, we have employed a simple spatial averaging in order to obtain homogeneous baryonic matter 
and have placed all instantons at the tip of the connected flavor branes, i.e., we have not determined their location 
in the bulk dynamically. All these deficiencies of our calculation can be improved in a systematic way, and thus we 
believe that our observations will help to move closer towards a holographic picture of real-world dense matter. For 
example, our observation of the first-order phase transition in the homogeneous approach is a clear indication that 
some of the physics away from the A = oo limit must be included, and we are currently working on extending the 
instanton gas approach in this direction. 

We have mainly focused on two crucial properties of dense matter, a first-order baryon onset and chiral restoration 
at large densities. Of course, fulfilling these necessary properties does not guarantee that we are quantitatively close 
to real-world QCD. Indeed, as we have argued, there are reasons to expect that, at best, the Sakai-Sugimoto model 
in the decompactified version can only be a rough guide to high-density QCD. Nevertheless, the lack of first-principle 
calculations and of any model for cold matter that is reasonably applicable over a wide density regime would make 
even a senri-realistic strong-coupling model very valuable. The Sakai-Sugimoto model in the limit considered here 
has only three free parameters, the ’t Hooft coupling, the Kaluza-Klein mass, and the asymptotic separation of the 
flavor branes. At this point, we have not attempted to fit them to any physical properties. After the model has been 
improved along the lines just indicated, it will be very interesting to fit them for instance to nuclear ground state 
properties and use the model for quantitative predictions, e.g., for the equation of state in the context of compact 
stars. 
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Appendix A: Confined geometry 

All derivations in the main text are done in the deconfined geometry. In this appendix we discuss the confined 
geometry which, in our context, mainly serves as a warm-up exercise for the numerical calculations of the deconfined 
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geometry. The reason is that it allows for a very simple scenario where the embedding of the flavor branes in the 
background geometry is trivial. Most of the arguments are completely analogous (in many instances simpler) as in 
the deconfined geometry. 

The induced metric on the flavor branes is 
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1. Baryonic matter from an instanton gas 


As for the deconfined geometry, we first expand the action for small non-abelian field strengths, use this expansion 
to introduce the instanton gas with the help of the instanton solution of the YM action, and then consider the action 
to all orders in the non-abelian field strengths. For the expansion we need 
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(A6) 

In order to discuss the instanton solution, we need to consider the YM action, which is immediately obtained from 
this expansion. We consider the case of maximally separated branes, such that the embedding of the flavor branes 
is trivial, x ' 4 = 0 (for a generalization to the non-antipodal case, see Ref. [75])• We also set the chemical potential 
to zero, i.e., ao(oo) = 0. Nevertheless, due to the CS term, do(u) acquires a nontrivial profile even in that case H2] , 
which affects the instanton solution. For now, we also ignore the CS contribution, such that we can set do (it) = 0. 
Then, dropping the term that is independent of the gauge fields, we obtain the YM action 
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where we have used Eqs. ([6]) and ( |A3| , the variable Z from Eq. ( fi~3| with U c replaced by Ukk (we now integrate over 
both halves of the connected flavor branes), and abbreviated k(Z) = t/]| K + Z 2 Ukk, h(Z) = (U^ K + ^ 2 ^kk) - 1 ^ 3 - 
In flat space, setting k(Z) = U^ K and h{Z) = L+ K . the YM equations of motion are solved by the instanton 
solution [42] 
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The reason for the appearance of this factor, which rescales Z relative to A', is the relative factor between the F~- and 
Ff z terms in the YM action (we have taken the flat-space limit by dropping the Z dependence, but not setting the 
constant prefactors to one). Since the instanton width will later be determined dynamically, we are free to denote it 
by p/ 7 . This will turn out to be convenient when we take the spatial average of the instanton gas, see Eq. ( |17[ ), after 
which 7 simply becomes a prefactor of the instanton profile. 

With F (jj, = d^A u — d v A^ + i[AA v ], the resulting field strengths are 
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The generalization of this single instanton to a gas of instantons is explained in the main text for the deconfined 
geometry, see Sec. Ill A In our approximation, this generalization results in the field strengths (21). Completely 
analogously, we obtain for the confined geometry 
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2u kk y/f mq{u) Tr [FijF ku ]eijk 


with 


3 * 4 K nj-f (2 na') 2 ’ ( M KK R ) 6 3m 1 / 2 (2 ira') 2 ’ ( M KK R ) 3 

9U 1 / 2 (p 2 wkk) 2 

q{u) = 


2 niq(u ) 
'3 {2na'y 


4 sfj ( u 3 - u^ K + p 2 mki<) 5/2 ’ 
With m kk defined in Table [I] and using Eq. (A3), we have 


M kk = g ' 


(A12) 


(A13) 


(A14) 


As in the deconfined case, our actual calculation is performed by keeping all orders in the field strengths. For the 
abelian case, this yields [compare with the corresponding expression of the deconfined geometry (29)] 


det(p + 27TOAF) = (Jj 


3/4 


(2W) 2 /4 + [1 + u 3 f 2 x' 2 + fa' 2 ] 

o , „ . s 1/2 

R\ ( 2 tt a'yfiFijFkudjk) 2 

U, 


, , R 

1 + t7 


(2 W) 2 F 2 


(A15) 


We now generalize this expression to the non-abelian case, following the prescription explained at the beginning of 
Sec. IIIB With the help of Eq. (A12), working with maximally separated flavor branes, x' 4 = 0, and including the 
CS term, we arrive at the Lagrangian 


C = 


u 5 ' 2 

VJ 


\/(l +ffi - /a , o 2 )( 1 + 92 ) - rua 0 {u)q(u), 


where 


,1/2 


Vf( u ) 


Vf( u ) 


( \ _ “ v j / \ ( \ _ “kkv j ( \ 

9iW = —^ 2 “^— niq(u ), g 2 {u) = —— rnq(u) 


with q(u ) from Eq. (A13). The solution of the equation of motion for do is given by 

( niQ ) 2 
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_ 1 + Pi 

W 5 / 1 + 02 + ^ 


(A 16 ) 
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where 
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Following the arguments in Sec. Ill B we minimize the free energy with respect to n/ and p, which leads to the 
coupled equations 


where 


0 = 


0 = 




du 


du 


— (^C 1 + —C ]+d'Q 

2Vf \dni dnj ' 0 


u 5 ' 2 

2V7 


(^ rl + ^ C ) +n/ “° 


dQ 

dp 


l l + gi-fa' 2 zt 5 / 2 Vfo 


\/l +9i 


1+52 


niQ 


1 + 52 + 


and we can write the free energy as 

n r°° , u 5 / 2 
77 


r°° u 5/2 _ / 

/ ^ u n \A + 5i \ 1 + 52 + 

JU KK V J ’ 


(n/Q) 2 


+/Q) 2 
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(A20a) 

(A20b) 

(A21) 

(A22) 


2. Baryonic matter from a homogeneous ansatz 


In the confined geometry with maximally separated flavor branes, the homogeneous ansatz for the gauge fields (56) 
yields, in analogy to Eq. (61), the action 


S = AT du 

J »KK 


.5/2 


where 
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\j(7+9i- /a'o 2 )(l +52) - ^Y~a 0 h 2 ti 

3Ag h(u) A 


V7 
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(A23) 


(A24) 


We now proceed completely analogously to the deconfined geometry, as explained in Sec. |IV A[ the baryon density is 
related to the boundary value of h in the bulk, nj = — 3Ao/2 /i(mkk) 3 , and the equation of motion for ao in integrated 
form is 


tt 5/2 v7aoVl +52 

V 1 + 51 - fa\J 


o\ 

= n I Q, Q(u) = 1 + — —h(u) 3 , 

2ni 


such that 


With the abbreviation 
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the equation of motion for h becomes 


d u 


5/2 C 1 dgA _ 9X 0 h 2 Cn I Q = u 5/2 C dg 2 
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(A25) 
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and we find that h behaves around u = ukk as given in Eq. (70) for the deconfined geometry. The minimization with 
respect to nr yields 
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and the free energy becomes 


fi 

A 7 


r°° u 5/2 _ / 

/ du-=-y / l + g 1 \l + g 2 + 

J w K K V / * 


(■niQy 


Hn /. 


(A30) 


[We have checked that minimizing this expression numerically with respect to nj gives the same result as using the 
simpler (A291.] 


Appendix B: Pointlike baryons, mesonic phase, and chirally restored phase 

In this appendix we discuss the baryonic phase in the approximation of pointlike baryons as well as the mesonic 
and chirally restored phases, which are needed for the discussion of the phase structure in the main text. Most of the 
discussion in this appendix concerns the deconfined geometry, but we also briefly discuss the (much simpler) results 
for the baryonic and mesonic phases of the confined geometry with maximally separated flavor branes. The chirally 
restored phase exists only in the deconfined geometry. 


1. Pointlike baryons 

a. Deconfined geometry 


As argued in the main text, we recover the Lagrangian for pointlike baryons (261 in the limit of vanishing baryon 
width p —>■ 0, 


C = C 0 + n I |\//t(w) - a Q (u) S(u-u c ), 


(Bl) 


with Co from Eq. (25). This Lagrangian was the starting point in Ref. |53| . and here we recapitulate the calculation 
of that reference in a concise way (a very similar calculation can be found in Ref. |55j . where a background magnetic 
held was included). The equations of motion in integrated form are 


u 5 / 2 a'n 


a /1 + vfifrx'l - dg 2 
u 5 / 2 u 3 f T x 4 

a/1 + u 3 f T x'l - d' 0 2 


= raj, 


= k, 


(B2a) 


(B2b) 


with an integration constant k that has to be determined in the following. These equations are easily solved for a' 0 
and x 4 , 


a' 0 (u) = 
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3/2 


\/ fr(i 
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x 4 (u) = 
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(B3) 


where we have abbreviated 


g(u) = (u 8 + ri 3 ?jj)/ T (ra) - k 2 . 


(B4) 


Consequently, do and x 4 are 


a 0 (u) = nj 



v 3/2 \Z7r(vj 


+ a 0 (u c ), 


x 4 (u) = k 



dv 

v 3/2 \JfT[v)\Jg{v) ’ 


(B5) 


where we have set x 4 (u c ) = 0. 

As explained in the main text for finite-size baryons, we need to minimize the free energy with respect to nj and 
u c . From the minimization with respect to nj we obtain 

do {uc) = y\//rW ■ 


(B6) 
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As a consequence, the on-shell action is solely determined by Co- The minimization with respect to u c yields [using 
Eq. (|B 6 ])] 


ri! 1 + /t(«c) _ Vdiuc) 


3/2 


This relation can be used to determine k as a function of u c and n/, 

k 2 = Oc + «c«/)/t(m c ) - u 3 c (y) 


( ni Y 

1 + /t(«c) 

V 3 ) 

2 


1 2 


(B7) 


(B 8 ) 


By inserting this result into the expression for x' 4 (u) from Eq. (B3), we see that x r 4 (u ) assumes a finite value at 
u = u c for all nonzero nj. Therefore, the embedding of the connected flavor branes acquires a cusp, in contrast to 
the instanton gas approximation, as shown in the main text. 

The quantities u c and nj are then determined by the coupled system of equations, 




du 


m 3 / 2 \/ f T {u)y/g{u) ’ 


r 

M - a 0 (u c ) =nj 

J U r 
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.3/2 


V fr{u) 


\fg(v) 


(B9) 


where we have imposed the boundary condition = ao(oo). Since the right-hand side of the second equation is larger 
than zero, there can only be solutions for p, > do(u c ). In the limit n/ —> 0, /. i assumes this minimum value, i.e., with 
Eq. (B 6 ), the critical chemical potential for the second-order onset of baryonic matter is 


Monset — g \/.[t(Uc) , 


where u c at nj —^ 0 is determined from 

K 


du ■ 


»c\//t(Wc) 


m 3 /2 \/u 8 fr{u) - u%f T {u c ) 
At zero temperature, this can be evaluated analytically, and we find 
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Monset (T — 0) 
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0.17495 

l 2 
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(BH) 


(B12) 


By solving Eq. (Bll) numerically we can compute the baryon onset for all temperatures. The result is the thick 
dashed line in Fig. [1J For the free energy comparison with the mesonic and chirally restored phases we need the free 
energy 
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pointlike 


= a r 


i 5 1 

auu - 


,3/2 




Vg(u) 


(B13) 


As mentioned in the main text below Eq. (39), a divergent vacuum contribution has to be subtracted, as for all free 
energies in this paper. With this simple renormalization and using the results from the following two subsections, 
one finds that the free energy (B13) is smaller than the free energy of the mesonic and chirally restored phases as 
soon as baryons are allowed to appear. Therefore, the chemical potential (BIO) marks the transition from vacuum to 
baryonic matter. 


b. Confined geometry 


The Lagrangian for pointlike baryons with maximally separated flavor branes in the confined geometry is 

~Ukk 


C = C 0 + m 


L 3 


- 5o(mkk) S(u — ukk) ■. 
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u 5 / 2 
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The equation of motion for do in integrated form is thus 

M 5 / 2 v7do 

V 1 - fa\? 


By defining the new variable 


y(u) = f 

J Uy 


= m . 


dv 


VfV v5 + n j ’ 


the equation of motion becomes 

d y d 0 = ni => a 0 (y) = niy + c. 

We determine the integration constant c with the help of the boundary condition do ( 2 / 00 ) = /T where 

du 


Voo = y{u = 00 ) = / 

J Uy 


\fjy/u 5 + n 2 j 


We find c = y, — niyoo and thus 

do (y) = n I {y-y oa ) +/i. 

Minimization with respect to n/ yields [completely analogously to the deconfined geometry, see Eq. (B6)] 
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As in the deconfined geometry, the free energy only receives a contribution from £o- With 
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(B23) 
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2. Mesonic Phase 


a. Deconfined geometry 


For the mesonic phase, we simply set nj = 0 in the Lagrangian (26 1 , such that C = Cq. This leads to the equations 
of motion in integrated form 


U 5 / 2 d'„ 


\/l + VpfTX'l - 


= 0 , 


U 5/2 U 3 frX 4 1 / 2 , , 4 

—===== = f T {u 0 )u 0 , 

V 1 + M 3 /t^ 4 - a 0 


(B25a) 


(B25b) 


where we have denoted the location of the tip of the connected branes by uq (at the second-order onset of pointlike 
baryons, uo = u c ), and where we have determined the integration constants on the right-hand side from the boundary 
conditions 


ao( u o)= 0 , d 0 (oo) = /i, ^ 4 (^ 0 )= 00 , 


du X 4 . 
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Consequently, 


Uq — 0 ; £4 — 


f 1/2 

Jt 


(v>o)l 


fT (u)u 3/2 y/u 8 f T {u) - Wo/t(wo) 


(B27) 


In particular, ao is constant, indicating that the baryon density vanishes in this phase. It thus remains to determine 
Mo as a function of temperature (it does not depend on //). This is done with the help of the last relation of Eq. (B26l. 
At zero temperature, where fr{u) = 1, one finds the analytical solution 


m 0 (T = 0) = ^ 
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The free energy becomes 
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Again, at T = 0, this can be evaluated analytically. Introducing a cutoff A, we obtain 


(B28) 


(B29) 


n 


mesomc 


(T= 0) 


r A 


U 


Af 


= / duu 5 / 2 - _ 

Ju 0 V U 8 - Mg 


2.r/2 _ 2 1 o 7t 4 T (‘ 16 ) tan 16 

7 15 F r (ff) 

""-V- 


Eli) 

r (ih) 


~ 5.4 x 10 ~ 3 /£ 


This value is used for example in the upper right panel of Fig. [ 6 ] 


(B30) 


b. Confined geometry 


Setting n/ = 0, the Lagrangian is simply given by £ = £0 with £0 
motion for oq in integrated form 




V 1 - fa 0 


(B31) 


Here we have set the integration constant to zero to ensure Mq(mkk) =0. As a consequence, the solution is trivial, 
do{u) = n, and the free energy becomes 


n = Af 


f°° , m 5 / 2 
I du ^j ■ 

UKK V J 


(B32) 


3. Chirally restored phase 


In the chirally restored phase, we start from the same Lagrangian as in the mesonic phase, but since the flavor 
branes are straight and disconnected, the holographic coordinate assumes values in the interval u £ [Mr,oo], and we 
have £4 = 0. Therefore, there is only one nontrivial equation of motion whose integration constant is the baryon 
density nj, 


u 5/2 a' 0 

V l ~a'c 


= ni. 


(B33) 


In this phase, there are no baryons, the baryon density is generated by quarks. With the boundary conditions 


ao(Mr) = 0, a 0 (oo) = //, 


(B34) 


we find 
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and the following implicit equation for nj as a function of fi and T, 


116 
5’ 2’ 5’ ~ rij 

At zero temperature, we have the simple result 
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The free energy is 
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By comparing this free energy to the one from the mesonic phase ( |B30| , we compute the critical chemical potential 
for the chiral phase transition at T = 0, 


Hc{T = 0 ) = 


167T 1 11 / 14 
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(B41) 


For nonzero temperatures, the critical chemical potential has to be computed numerically. The result is shown in 
Fig- 0 where we see that for small temperatures there is no chiral phase transition because the mesonic phase is 
superseded by the baryonic phase before the chirally restored phase becomes favored over the mesonic phase. 
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